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Overlap Fermions on a 20 4 Lattice 
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We report results on hadron masses, fitting of the quenched chiral log, and quark masses from Neuberger's 
overlap fermion on a quenched 20 4 lattice with lattice spacing a = 0.15 fm. We used the improved gauge action 
which is shown to lower the density of small eigenvalues for H 2 as compared to the Wilson gauge action. This 
makes the calculation feasible on 64 nodes of CRAY-T3E. Also presented is the pion mass on a small volume 
(6 3 x 12 with a Wilson gauge action at (3 = 5.7). We find that for configurations that the topological charge 
Q / 0, the pion mass tends to a constant and for configurations with trivial topology, it approaches zero possibly 
linearly with the quark mass. 



Overlap fermion has great promise in studying 
chiral symmetry on the lattice |l]] . Recent numer- 
ical implementation of the overlap fermions has 
led to studies of chiral condensate quark 
mass {I, as well as the checking of chiral sym- 
metry and scaling (JH. However these studies 
are limited to small volumes due to the enor- 
mous numerical cost associated with approximat- 
ing the matrix sign function which is about fifty 
times more than inverting the quark matrix of the 
vanilla Wilson action. 

In the optimal rational approximation Q of the 
matrix sign function of Neuberger's overlap oper- 
ator @ 
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where e(H) = H/ V H 2 and H is taken to be the 
hermitian Wilson-Dirac operator, i.e. H = 75-Dm 
with 0.25 > k > k c ), it is cost effective to project 
out a relatively few eigenmodes with very small 
eigenvalues in the operator H 2 in order to reduce 
the condition number and speed up the conver- 
gence in the inner do loop MM . At the same time, 
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it improves the chiral symmetry relation such as 
the Gell-Mann-Oakes-Renners relation Q. How- 
ever, it is shown ||] that the density of these small 
eigenmodes grows as with a being the lattice 
spacing. As a result, it is very costly and imprac- 
tical to work on large volumes with the currently 
used lattice spacings. There are simply too many 
small eigenmodes to be projected out. 

For this reason, we explore other options to 
clear this hurdle. We have tested the 0(a 2 ) im- 
proved gauge action of Morningstar and Pear- 
don H and find that the density of these small 
eigenvalue modes is decreased to a point that it 
becomes feasible to go to large volumes with a size 
of the lattice 3 to 4 times of the Compton wave- 
length of the lightest pion. We further find that 
the anisotropic action |J requires projection of 
more small eigenvalues in H 2 in order to achieve 
the same convergence in the inner loop than does 
the isotropic one. Thus, we use the isotropic ac- 
tion. We also find that using the clover action 
with either sign requires the projection of more 
small eigenvalue modes. Therefore we use the 
Wilson action for H in the Neuberger operator 
with k = 0.19. On a 20 4 lattice with a = 0.15 
fm as determined from the Sommer scale rn, we 
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project out 80 small eigenvalues. Beyond these 
eigenmodes, the level density becomes dense. As 
a result, the number of conjugate gradient steps 
for the inner loop is about 160 and the conju- 
gate gradient steps is about 210 for the outer 
loop. These numbers are about the same as those 
for the Wilson gauge action on small volumes ||] . 
Therefore, other than the overhead of projecting 
out the small eigenvalue modes, the cost scales 
linearly with volume. 

We shall report preliminary results on 38 gauge 
configurations. Given a = 0.15 fm for the 20 4 
lattice, the physical length of the lattice is 3 fm. 
We shall study pion mass as low as ~ 240 MeV 
so that the size of the lattice is about 3.6 times 
of the pion Compton wavelength. 

We first show the pseudoscalar mass squared 
m 2 P a 2 as a function of moa in Figure We fit 
them in the following form jl0| with and without 
the quenched chiral log term 6 

m 2 P a 2 = c+2Am a 2 {l-Sln(2Am /A 2 x )}+4:Bml.(2) 
20 4 improved gauge, overlap 
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Figure 1. Pion mass squared as a function of the 
bare quark mass moa on the 20 lattice with a = 
0.15 fm. 



As shown in Table [|, the fits with and without 
the quenched chiral log S give comparable x 2 /DF 
for the 6 to 10 smallest quark masses ranging from 
0.014 to 0.24 (see Figure pi). The fit is insensitive 



to A x in the range of 0.6 to 1.4 GeV. We list 
results with A x = 1.0 GeV in Table @. Further- 
more, the errors on S are much larger than the 
corresponding central values. We conclude that 
there is no signal for the quenched chiral log in the 
range m^/mp ~ 0.35 — 0.80. This is in contrast to 
the study of the Wilson fermion action on lattices 
with the same physical volume and comparable 
quark masses Q which supports the presence 
of quenched chiral logarithms with a magnitude 
5 ~ 0.1. We also note that in our fit, the intercept 
c is consistent with zero which is to be expected. 

Next we plot the nucleon, the vector and pseu- 
doscalar masses as a function of mo a in Figure |2|. 
A simple linear fit with 10 smallest quark masses 
gives mya = 0.534(12) at the chiral limit with 
X 2 /DF = 0.05. Using the r to set the scale, this 
corresponds to 712(16) MeV. Since the error in 
tojv is still large, we can not draw any conclusion 
on the misr/mp ratio yet. 
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Figure 2. Nucleon, vector and pseudoscalar me- 
son masses. 



We have computed the quark mass through 
the chiral Ward identity ZaQ^A^ = 2Zg 1 moZpP 
where = ^7^75(1-/2)^ and P = 4ji^^{t /2)ip. 
We plot |£gJmp( M ) = Z^Mmo in Figure | 
and fit them in the form c + Am§a + Bm\a 2 . The 
fits with a few smallest quark masses are given 
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Table 1 



The fitted parameters in Eq. (|2 


|). (5 = corresponds to the case without the chiral log. 




#pts. 


c 


S 


A 


B 


X 2 /DF 


5 


0.0016(50) 





4.51(22) 





0.008 


6 


-0.0006(41) 





4.65(13) 





0.2 


6 
6 


0.0045(70) 





4.17(56) 


12(24) 


0.02 


8 


0.0047(53) 





4.15(30) 


7.5(29) 


0.03 


8 


0.004(18) 


0.015(297) 


4.13(56) 


8.2(149) 


0.035 


10 


0.0051(49) 





4.12(22) 


7.8(20) 


0.04 


10 


0.02(15) 


0.043(227) 


4.07(36) 


9.7(98) 


0.04 



Table 2 

The renormalized quark mass §jf~j m^ s (fj>)a fitted in the form c + Arriba + Bm§a 2 . 



# 


c 


A 


B 


X 2 /DF 


Z A (ji,0) 


4 


-0.00041(34) 


0.939(24) 


1.66(34) 


0.003 


1.065(27) 


5 


-0.00049(24) 


0.945(14) 


1.56(16) 


0.05 


1.058(16) 


6 


-0.00063(16) 


0.956(8) 


1.43(7) 


0.29 


1.046(9) 



in Table ||. We see that the intercept c is consis- 
tent with zero for the case with 4 smallest masses. 
From this we can deduce the non-perturbative 
renormalization constant Za{h, mo) as a function 
of moa, i.e. (p,,mo) = A + Bm^a. We see 
from Table |^ that Za{^, tuq = 0) is fairly far from 
the tree-level value of 1.368. 

Finally we shall explore the behavior of the 
pseudoscalar meson mass on a small volume. 
Starting from the generalized Gell-Mann-Oakes- 
Renners relation 

m a J d 4 a;(7r(x)7r(0)) = 2(V^), (3) 

and assuming pion dominance, one obtains the 
usual Gell-Mann-Oakes-Renners relation 

ml= 2 ^m. (4) 

In the quenched approximation, the chiral con- 
densate has the following behavior for small mo 

(#>^ + a + &m , (5) 

where the first term due to the fermion zero 
modes associated with the topological charge Q 
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Figure 3. The quark mass in the MS scheme. 
The solid line is a linear plus quadratic fit with 
the 4 lightest mass points. 



is the quenched artifact which is prominent as 
mo — > and/or at small volume V. This term 
is observed in {ipip) through a direct calculation 
with the overlap fermion Thus, barring ad- 
ditional complication due to the quenched chiral 
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6 3 x 12 with f3 = 5.70, 50 configurations 
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6 3 x 12 with /3 = 5.70, Q = 0, 18 configurations 
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Figure 4. m,p S a 2 as a function of mo a. 



log, the pion mass is expected to have the follow- 
ing dependence on the quark mass mo 
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We display in Figure f| m\a} as a function of 
mo a for the 6 3 x 12 lattice with 50 Wilson gauge 
configurations at /3 = 5.7. The top one shows that 
m\ approaches a constant in the chiral limit as 
shown in Eq. (^J) . To further verify that the con- 
stant is indeed due to the nonvanishing topologi- 
cal charge in the configuration, we separate the 50 
gauge configurations into 32 Q ^ ones (the mid- 
dle figure) and 18 with Q = (the bottom figure). 



for the Q ^ case approaches 
the same constant as mo<z — ► 0, while it is much 
smaller for the Q = case which has a tendency 
to approach zero at the chiral limit. From the 
curvature of the dependence on m a, one spec- 
ulates that moa is below the Thouless energy in 
such a small volume so that the amo term in Eq. 
(||) vanishes. As a result, m\ is proportional to 
m,Q and consequently the pion mass approaches 
the chiral limit linearly with mo. However, one 
needs more statistics to confirm this scenario. 

To conclude, we demonstrate that it is feasible 
to implement the overlap fermions on large phys- 
ical volumes. We do not see the quenched chiral 
log in the range m^/mp ~ 0.35 — 0.80. 

This work is partially supported by U.S. 
DOE Grants DE-FG05-84ER40154 and DE- 
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